Fractional orbital occupation of spin and charge in artificial atoms 
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We present results on spin and charge correlations in two-dimensional quantum dots as a function 
of increasing Coulomb strength (dielectric constant). We look specifically at the orbital occupation 
of both spin and charge. We find that charge and spin evolve separately, especially at low Coulomb 
strength. For the charge, we find that a hole develops in the core orbitals at strong Coulomb 
repulsion, invalidating the common segregation of confined electrons into an inert core and active 
valence electrons. For excitations, we find a total spin-projection Sz = —1/2 breaks apart into 
separate occupations of positive and negative spin. This dissociation is caused by spin correlations 
alone. Quantum fluctuations arising from long-range Coulomb repulsion destroy the spin dissociation 
and eventually results in all orbitals carrying a negative spin. 
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PACS numbers: 85.35.Be, 73. 21. La, 75.75.-|-a, 73.21.-b 



Quantum dots, especially at zero magnetic field, can 
often be considered artificial atoms [ij, Like real 
atoms, artificial atoms can confine a definite number 
of electrons, and the electronic and magnetic proper- 
ties are largely governed by the actual number of con- 
fined particles. On the other hand, there are dramatic 
differences between the two which stem primarily and 
rather obviously from the synthetic nature of quantum 
dots. Quantum dots are not identical in the quantum- 
mechanical sense, and they are more tunable 3] than 
most atomic and condensed matter systems. Both the 
qualitative and quantitative features of the confinement, 
for example, are tunable; this, in turn, can be exploited to 
control the angular- momentum 4] and the spin 0, 0, 0] 
content of ground and excited states. With respect to 
the present work, the most significant difference between 
atoms and electrostatically-defined quantum dots in par- 
ticular is their size. Due to their much larger size, 
quantum dots typically have single-particle level spacings 
of order 1 meV — approximately 10"* times smaller than 
atoms. Consequently, the relative strength of Coulomb 
repulsion is also 10"* times larger in quantum dots than 
it is in atoms. In these quantum dot systems. Coulomb 
effects are never negligibly small and, for few electron 
systems especially, the long-range nature of Coulomb re- 
pulsion often drives the observed phenomena 

Our focus in this work is on the orbital occupation of 
particles. We seek to determine which orbitals are occu- 
pied by charge, which by spin, and how much of each is 
contained in the given orbitals. We have investigated this 
question as a function of decreasing dielectric constant, 
effectively increasing the strength of long-range Coulomb 
repulsion. We generally find that correlations can be in- 
duced both by the spin symmetry as well as Coulomb ef- 
fects. Spin correlations are most prominent at relatively 
weak interactions and we find in this regime marked dif- 
ferences between the distribution of charge on the one 
hand and spin on the other. Coulomb fluctuations gen- 



erally destroy this spin-charge distinction, resulting in 
identical distributions of both spin and charge. 

We specifically report here on three main results. First, 
the common segregation of the charge into an inert core 
filling the low-lying shells and an active valence shell 
completely breaks down at strong interactions; the core 
shells are neither filled nor inert. We find a hole de- 
veloping in the core as electrons maximise correlations 
by occupying higher shells, depleting the core. Second, 
we find that the evolution of the fractional occupation 
of spin is different than that of charge, making it useful 
to consider spin and charge as independent entities even 
though they are both attached to each electron. Thirdly, 
we look at excitations which preserve the symmetries of 
the ground state and find significant correlation effects 
due entirely to the spin symmetry, even in the complete 
absence of Coulomb effects. We find that a spin excita- 
tion of Sz — —h/2 does not merely delocalize across the 
various orbitals, but rather dissociates into a relatively 
large negative spin and a small positive spin, such that 
the global spin symmetry is preserved. Quantum fluctu- 
ations due to long-range repulsion destroy this effect and 
lead to a single delocalized spin with Sz = ~h/2. 

Our model is a simple two-dimensional system of N 
electrons confined to a circular parabolic potential with 
long-range Coulomb repulsion. The Hamiltonian is 
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where the dielectric constant e = ecaAs/a is scaled by a 
relative to that of GaAs (ecaAs ~ 12.4). Thus, a = de- 
scribes the non-interacting limit with Coulomb repulsion 
growing linearly with a, and with a = 1 describing the 
GaAs system. Our results below are expressed as func- 
tions of a. In all cases with a ^ 0, we consider the full 
long-range (~ 1/r) repulsion. 



2 



Our strategy is to first obtain the many-body eigen- 
states of Eq. (Q), including correlation effects. We sub- 
sequently analyze these many-body states (rather than 
just the energies) and directly determine the orbital oc- 
cupations of the low- lying states. 

The numerical approach we adopt is guided by three 
main requirements. First, we require the full many-body 
states, properly accounting for all exchange and correla- 
tion effects; techniques based on effective single-particle 
orbitals or those based on obtaining energies only are 
therefore insufficient. Second, we also require excitations 
above the the ground state. Indeed we show below that 
an excitation of a single Sz ~ ^1/2 system results in 
both positive and negative spin occupation of orbitals. 
Thirdly, in these confined systems, the strongest corre- 
lations often occur for small occupation numbers where 
appeals to screening effects are invalid. 

For these reasons the most appropriate numerical tech- 
nique is a configuration-interaction (exact diagonaliza- 
tion) approach. This affords the full many-body corre- 
lations of the ground state and the excitations, and the 
few-body systems we consider are not beyond the reach 
of techniques with a heavy computational load. 

The Hamiltonian Q is exactl y so lvable in the non- 
interacting (a = 0) limit 0, lul llll|. In this case, the 
single-particle eigenstates are a pair of harmonic oscilla- 
tors |nm), whose notation is given in the quantum 
number n is a Landau level index, whereas m denotes a 
guiding center. The angular momentum of orbital \nm) 
is Lz — {n — m)h and the energy at zero field of the same 
orbital is i?„i„ = fkuoin + m + 1). 

Including spin, then, our single-particle orbitals are 
denoted by Inms). We build the many-body ba- 
sis states from antisymmetrised tensor products of 
the single-particle eigenstates, jnirriiSi, n2m2S2, • ■ ■) = 

The various rotational symmetries of Eq. |^ imply 
conservation of angular momentum projection L^, spin 
projection Sz, and spin magnitude S^. We use these 
quantum numbers to label the correlated eigenstates of 
the interacting problem. Lz and Sz do not induce corre- 
lations, but S^ may induce correlations. Indeed we show 
this to be the case particularly for the excitations even in 
the absence of Coulomb repulsion. In this regime, spin- 
correlations are strongest. 

We denote by \Sz,Lz,i) the i"^ state 
jnimisi, n2TO2S2, • ■ ■} with fixed angular momen- 
tum Lz = '^kfi{nk — Tnk) and fixed spin projection 
Sz = J2k^^k/2 {sk = ±1). This state is a single (un- 
correlated) antisymmetrised state and is an eigenstate 
of Eq. m in the absence of interactions. Similarly, by 
\S,Sz,Lz,j), we denote the j'^ spin eigenstate with 
spin magnitude h'^S{S + 1). This state is in general a 
correlated state 

\S,Sz,Lz,j)=Y.al\Sz,Lz,i), (2) 
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whose correlations are due entirely to the spin symmetry. 
This correlated state is also an eigenstate of Eq. (Q) — for 
a = — and the spectra will generally consist of degener- 
ate subspaces. Finally, the state \E, S, Sz,Lz, k) denotes 
the /c*'' eigenstate of Eq. with energy and with 
arbitrary a. This state contains Coulomb correlations in 
addition to those imposed by the spin symmetry. It can 
be written as a coherent superposition of the (correlated) 
spin eigenstates in Eq. (O: 

|i?,^,^„L„A:) =^/3||5,^„L„j). (3) 

These states will in general lift the degeneracies contained 
in the spin-correlated states alone. 

Our immediate objective is to find all E and all (3^ in 
Eq. JSJ. In the configuration- interaction method, Eq. Q 
can be viewed either as a variational ansatz with the j3^ 
and the E as variational parameters, or, equivalently, as 
a direct matrix diagonalization in the basis of the states 
IS*, Sz,Lz,j), where the (3^ and the E emerge as the eigen- 
vectors and eigenvalues of the diagonalization procedure. 
In cither case, both views yield identical results and both 
yield the exact solution in the limit of an infinite number 
of terms in the sum of Eq. 

To construct our many-particle basis states, we start 
with a sufficiently large set of single-particle orbitals 
\nms). In these calculations, we have taken 288 single- 
particle states. For a three-particle system with Sz = 
— 1/2, this gives over 10^ possible three-body states. 
From these, we consider only those states, Eq. con- 
sistent with the desired quantum numbers. These sym- 
metry requirements dramatically reduce the number of 
basis states and, in practice, wc find it rarely neces- 
sary to consider greater than ^ 10"^ many-body states 
\Sz, Lz,i) ■ From these states, we form all correlated spin 
eigenstates \S, Sz, Lz, j), Eq. (0), for a given quantum 
number S. The can be determined by direct diago- 
nalization of the spin operator S'^ in the \Sz,Lz,i) ba- 
sis iX^. Alternatively, a closed, analytic, and exact form 
for the al can be determined given the number of singly- 
and doubly-occupied orbitals using group t heory meth- 
ods and Clebsch-Gordon technology [l3l [iJI . 

We use these computed correlated spin states, Eq. 
as a basis to directly diagonalise the Hamiltonian Q, 
thus obtaining the coefficients (3j in Eq. (jSJ as well as 
(less importantly in the present case) the energies E. A 
closed form expression for arbitrary many-body matrix 
elements of the long-range Coulomb interaction exists di- 
rectly in spectral space 0,^3]' obviating the need for real 
or Fourier space integration. 

For a given state jV') of the system, the charge occupa- 
tion (number of electrons) on the single-particle orbital 
|nm) is given by 

pT = (V'KcLtC^T + cLiC«mi)l^), (4a) 
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FIG. 1: Fractional orbital occupation of charge in the ground 
state. The color-scale plot shows how the three particles dis- 
tribute themselves across the various shells (n+m) at different 
interaction strengths. With increasing interaction strength, a 
hole develops in the core orbital which, at zero interaction, is 
doubly occupied. The plot depicts the three-particle ground 
state at zero field, S = 1/2, = -1/2, = -1. 

whereas the spin (Sz) occupation of the same orbital is 

Sr = (^KcLtCn^T - ci^iCnmim. (4b) 

Once the eigenstates have been computed, the expres- 
sions in Eq. are evaluated for each orbital (m, n). 

A three-particle system, at i? = 0, is perhaps the sim- 
plest system containing non-trivial spin and charge cor- 
relations. Since an odd-number of electrons has finite 
spin for all states, we can determine how both spin and 
charge distribute themselves. For weak interactions, the 
third electron is the sole occupant of the second shell; 
as Coulomb strength increases, charge correlations can 
be expected to play an increasingly significant role in 
the ground state and elementary excitations. A three- 
particle system, even at fixed orbital occupations |3 
may form correlations entirely due to the spin- physics, 
independent of Coulomb correlations. 

At weak fields, the ground state is formed from the 
manifold of states with definite total spin S = 1/2, spin 
projection Sz — —1/2, and angular momentum Lz = — 1 
(all in units of h) . 

For the charge degrees of freedom, we have the general 
sum rule X^mn P^j" — where iV = 3 is the number of 
confined electrons. At zero field, the single-particle spec- 
trum consists of degenerate shells, where the fcth shell 
{k = 1,2,3,...) consists of all orbitals \nm) satisfying 
n + m — fc — 1 . Figure ^ shows how the particles smear 
across the shells for various interaction strengths. At 
zero interactions, the exact ground state is a single anti- 
symmetrised state of a doubly occupied core (first shell) 
and a single electron on the second shell ■ As a func- 
tion of increasing interaction strength, we find that the 
distribution broadens. This is expected since the orbital 
states are no longer eigenstates of the Hamiltonian. Sur- 
prisingly, however, the peak of the distribution does not 
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FIG. 2: Separate spin and charge evolution with increasing 
interaction strength. The panel at left displays the change in 
the mean shell occupied by spin and charge separately. At 
right is shown the average angular momentum for both spin 
and charge. At zero interactions, the spin distribution signif- 
icantly differs from that of the charge, both among shells and 
angular momentum. Coulomb-induced fluctuations destroy 
this order and restore congruence among charge and spin. 



remain on the core shell, nor does the occupation of the 
core drop to unity, thereby allowing exchange to dom- 
inate. Instead the core occupation drops below unity, 
almost to zero, thereby leaving a hole in the core. The 
electrons gain both exchange and correlation energy by 
occupying higher shells at the expense of single-particle 
energy. 

In addition to the charge occupation, we can look to 
the spin occupation. The ground state we are considering 
here has = 1/2 and Sz = —1/2. At zero interactions, 
the doubly occupied core shell carries no spin; the spin 
sits in the second shell. Similarly to the charge, this ex- 
cess spin delocalizes across the orbitals with increasing in- 
teraction strength. The evolution of the spin, however, is 
distinct from that of the charge. In Fig. [3 the distinction 
is evident. The left panel displays the mean shell occupa- 
tion for the ground state charge, J2m ri("^ + ")Po'"' 
the ground state spin. Also shown, in the right panel, is 
the average angular momentum for both spin and charge. 
For the charge, the total angular momentum is set by the 
global quantum numbers to Lz — —h and so the average 
is fixed for all values of a. These symmetry constraints 
do not apply to the (excess) spin. Instead, the spin obeys 
the sum rule X)m.n = -^2. 

The distinction between spin and charge distributions 
is most evident at weak interactions. Here, only spin 
correlations play a prominent role. Upon increasing 
Coulomb strength, these additional fluctuations destroy 
this distinction, resulting in identical mean shell and 
angular-momentum occupations. 

We turn now to excitations above the ground state. 
In particular, we look at excitations which preserve both 
angular momentum and spin. We find even stronger cor- 
relation effects due entirely to the spin symmetry, inde- 
pendent of the Coulomb correlations. To see this, we look 
first at the noninteracting limit at zero field, and look at 
the lowest excitations with the same quantum numbers 
as the ground state, {N,S,Sz,Lz) = (3,1/2,-1/2,-1). 
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FIG. 3; Spin and angular momentum conserving excitations 
at zero interactions. At centre is the ground state. The top 
row depicts one-body excitations and the bottom row two- 
body excitations. The left column denotes spin correlated 
states (two per image) while the right denotes spin uncorre- 
lated states. 
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FIG. 4: Spin dissociation at zero interactions and its restora- 
tion with Coulomb repulsion. In the presence of only spin 
correlations {a — 0), the spin 5z — —1/2 breaks apart into 
a small positive component and a larger negative component, 
each occupying distinct regions of the quantum dot. Coulomb 
correlations destroy this dissociation, resulting eventually in 
a single delocalized Sz — —1/2 particle. The color bar is in 
units of ft/2. 



The lowest-energy excitations lie 2hLUQ above the 
ground state, and the subspace is six-fold degenerate. 
These excitations can be classified by two orthogonal 
criteria — whether they are a one-body or a two-body ex- 
citation, and whether they are a spin correlated or un- 
correlated state. The orbital configurations of these de- 
generate excitations are depicted in Fig. |21 The centre 
image depicts the ground state (uncorrelated) configura- 
tion. The top set of plots shows the two one-body excita- 
tions, formed, in this case, by moving a single particle up 
two shells. (The shell spacing is huJo-) The bottom set 
depicts equally energetic configurations, involving an ex- 
citations of two particles, each of which moves up a single 
shell. In contrast, the excitations can also be classified ac- 
cording to their correlations. The two images on the right 
are uncorrelated states; after antisymmetrization, these 
states are proper spin doublets {S = 1/2). There are no 
additional degrees of freedom due to spin, the states on 
the left, however, are not proper spin states even after an- 
tisymmetrization. These configurations, containing only 
singly-occupied orbitals, should be viewed as each repre- 
senting three states which dilfer solely by a permutation 
of the spins. From these three states, two orthogonal 
5 = 1/2 states can be constructed. (The third state 
has S = 3/2.) These doublet states are correlated; they 
cannot be written as a single antisymmetrised (Slater de- 
terminant) state in any basis. The correlations, however, 
are entirely induced by the spin-symmetry requirements 
and not from Coulomb interactions. (Indeed there are no 
interactions at a = 0.) 

The implication of correlated excitations is that the 
excess spin {Sz = —1/2) does not sit on a single orbital 
as in the ground state, even when there is no interaction 
among the electrons. In fact, it is not even the case that 



the excess negative spin merely delocalizes. This single 
negative spin dissociates into a small positive spin and 
a larger negative spin and these are each distributed on 
different orbitals. This dissociation can be traced back 
directly to the correlated doublets in the degenerate man- 
ifold. This is shown in Fig. 21 where we plot the net 
spin (5*2) per shell as a function of increasing Coulomb 
strength. A region of positive spin exists in the third 
shell while all others have a net negative spin. (The total 
spin is always Sz — —^2.) We stress that the dissocia- 
tion of the spin is not due to Coulomb correlations but 
is entirely due to the spin correlations. In fact, turning 
on the Coulomb correlations destroys the effect. This 
is also shown in Fig. ^ the quantum fiuctuations asso- 
ciated with long-range Coulomb repulsion act to destroy 
the positive spin contribution until all the remaining spin 
is only of one species. 

Correlation effects can be due to either long-range 
Coulomb repulsion or spin symmetry. At low fields, 
where low spin is often energetically favoured, spin- 
correlation effects are strongest and novel spin-induced 
phenomena may occur. Coulomb fluctuations in general 
induce different, often much stronger, effects and these 
charge correlations often overwhelm, destroy, or other- 
wise mask effects due to the spin symmetry. We have 
shown here that in the relatively simple and clean system 
of three interacting electrons with both orbital and spin 
rotational symmetry the distinction and dissolution of 
spin-symmetry-induced and Coulomb-induced phenom- 
ena can be demonstrated at moderate Coulomb strength. 

This work was supported by NSERC of Canada and 
by the Canadian Foundation of Innovation. 
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